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ABSTRACT 

Wc investigate the emission properties that a large-scale helical magnetic field imprints 
on AGN jet synchrotron radiation. A cylindrically symmetric relativistic jet and large- 
scale helical magnetic field produce significant asymmetrical features in transverse 
profiles of fractional linear polarization, intensity, Faraday rotation, and spectral index. 
The asymmetrical features of these transverse profiles correlate with one another in 
ways specified by the handedness of the helical field, the jet viewing angle {9ob), and the 
bulk Lorentz factor of the flow (F). Thus, these correlations may be used to determine 
the structure of the magnetic field in the jet. In the case of radio galaxies [6ob ^ 1 
radian) and a subclass of blazars with particularly small viewing angles (6*06 ^ 1 /F) , 
we find an edge-brightened intensity profile that is similar to that observed in the 
radio galaxy M87. We present observations of the AGNs 3C 78 and NRAO 140 that 
display the type of transverse asymmetries that may be produced by large-scale helical 
magnetic fields. 

Key words: galaxies: active, galaxies: jets, (magnetohydrodynamics) MHD, radiation 
mechanisms: non-thermal, polarization 



1 INTRODUCTION 

Large-scale magnetic fields are thought to play a central 
role in launching relativistic jets in active galactic nuclei, or 



AGN (|Blandford fc Znajek|[l977| [Blandford fc Payne|p82 



[Lovelace et al. 19871. Following the initial production of 
the jets, the large-scale magnetic field naturally develops 
into a helical configuration in which the toroidal compo- 
nent of the field coUimates the outflow into a narrow jet 
via magnetic hoop stresses (Benford 1978 Chan & Hen- 
riksen|[T980 1 . While the importance of large-scale magnetic 



fields in the initial production of jets is rarely disputed, it 
is not known how far from the nucleus the ordered large- 
scale component of the field continues to remain dominant 
over the tangled component. Current-driven instabilities, of 
which the m = 1 kink mode is most important, are prime 
candidates for tangling or mixing the large-scale helical field 



when it is toroidally dominated (Begel man||1998| jCiannios 
fc Spruit||2006 1. The kink mode instability may imply that 



the field becomes tangled soon after the large-scale field has 



launched, accelerated, and coUimated the jet (e.g., Marscher 
et al.'2008) , or it may imply that only tangled fields ( Heinz & 



Begelman„2000) or only poloidal fields ( |Spruit et al.||1997 l 



play a significant role in astrophysical jets. Alternatively, 
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the growth of the kink mode may be stabilized through 
various effects such as gradual shear or an external wind 



( Hardee||2004| [McKinney fc Blandford|[2009| and references 
therein). Clearly, observations are needed to further illumi- 
nate whether parsec scale AGN jets contain large-scale he- 
lical magnetic fields. 

The most direct method of observing the existence and 
geometry of large-scale magnetic fields, very long baseline 
interferometry (VLBI) polarimetry, yields largely ambigu- 
ous results. Parsec scale radio jet synchrotron emission is 
often highly polarized, implying that the jet magnetic field 
is anisotropic but not necessarily large-scale (Pacholczyk 
|1970[ ). The electric vector position angle (EVPA) of a syn- 
chrotron emitting element is perpendicular to the element's 
magnetic field direction projected onto the sky. However, 
due to relativistic aberration in AGN jets the EVPA is not 
necessarily perpendicular to the observer frame jet magnetic 



field (Lyutikov et al. 20051. Statistically, observed EVPAs 



have a bimodal distribution: most are either aligned with 



the local jet axis or perpendicular to it (Cawthorne et al. 
19931 'Marsch er et al.|[2502l [Lister fc Homan||2005[ [Kharb 
et al.||2008J . Two classes of models explain these data: (i) 
models assuming a helical field and a cylindrical jet (Lyu- 
tikov et al.||2005" l and (ii) shock and velocity shear models 
which assume that shock compression along the jet axis pro- 
duces EVPAs parallel to the jet axis and that velocity shear 
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explains EVPAs perpendicular to the jet axis (Laing 
Hughes et al.||1989a|b| [T99T| [Attridge et lal|1999| |Kharbl 
et al.|2005| )." Both models explain the presence of linear po- 



larization, which would be absent if the magnetic field were 
tangled and isotropic. 

To distinguish between these two explanations of VLBI 
polarimetry data and more generally test the assumption of 
parsec scale helical magnetic fields, several researchers have 
calculated how intensity and fractional polarization change 
across a jet in the transverse direction (i.e. profiles). Laing 



et al. ( 1981 2006 1 compared the polarization and intensity 
profiles of kpc jets to a variety of helical field models and tan- 
gled anisotropic field models and found the predicted asym- 
metries of helical models to be inconsistent with observed 
symmetrical profiles. However, [Papageorgiou ] p006| ) found 
a number of similarities between the intensity and polariza- 
tion profiles due to helical fields calc ulated in^Laing (1981 ) 
and AGN jet observations. Similarly, |Aloy et al.| ( |2000[ ) car- 
ried out three-dimensional hydrodynamnic jet simulations 
that found asymmetries similar to those discussed in this pa- 
per in profiles of intensity and polarization. |Lyutikov et al.| 
( 2005 1 studied the polarization profiles and other polariza- 



tion properties of relativistic jets with helical fields, finding 
that the presence of helical fields explains both the bimodal 
distribution of EVPAs for unresolved jets and the abrupt 
90° flips in the EVPA for jets with small viewing angles (i.e. 
blazars) . 

Recently, observed gradients in Faraday rotation mea- 
sure (RM) have received much attention. Following the sug- 
gestion that helical fields cause RM gradients across jets 
( Laing|[l981 Blandfordl|1993 1, observers have firmly estab- 



lished such behavior in the famous quasar 3C 273 ( |Asada| 



dients in other parsec scale AGN jets ( 


Asada et al. 


2008b 


Gabuzda et al.|2004||2008||Kharb et al.||2009, ,0'Sul 


ivan &| 


Gabuzda||2009 |Croke et al.||2010 


1. These observations have 



prompted a number of efforts to fit the observed RM pro- 
files to simple analytic models of AGN ( [Contopoulos et al.| 
2009 Konigl|2010 1 as well as numerical simulations of AGN 



( Broderick fc McKinney|2010 1 



To investigate the question of large-scale helical fields, 
we calculate transverse profiles using the polarized syn- 
chrotron absorption and emission coefficients, retaining the 
coefficients' dependence on the angle between the jet frame 
line of sight and magnetic field. These calculated profiles 
show that a jet having an axially symmetric structure and 
helical magnetic field nevertheless display asymmetric pro- 
files in intensity, Faraday rotation, fractional polarization, 
and spectral index due to the large-scale helical structure of 
the magnetic field. In this work we do not discuss the ob- 
servations or modeling of EVPA swings of bright emission 
features in parsec-scale AGN jets which may help further 
our understanding of AGN jet magnetic fields (Marscher et 
ar||2008| [Abdo et al.||2010[ ). This paper is structured as fol- 
lows. In ^ we discuss the assumptions made regarding the 
structure and geometry of the relativistic jet. Optically thin 
intensity profiles are calculated and compared to observa- 
tions in !j3] Parsec scale RM calculations and observations 
in ^ are followed by unconvolved and convolved calcula- 
tions of polarization profiles and polarization observations 
in 35| Spectral index profiles are calculated in !|6] Finally, 
in 97]we summarize our work and discuss how asymmetrical 



features in each of the aforementioned profiles should cor- 
relate with one another. As an example of how to employ 
the predicted correlations of our helical model, in S|7]we also 
briefly discuss possible correlations seen amongst the differ- 
ent transverse profiles of the radio jets 3C 78 and NRAO 
140. 



2 A SIMPLE MODEL 

2.1 Magnetic Field Structure 

A crucial assumption we make regarding the structure of 
the helical magnetic field is that the jet frame azimuthal 
and axial fields are of comparable magnitude, B'^/B'^ ~ 1 
(all primed quantities refer to the rest frame of the jet). 
As explained in §2.2| if this assumption is false and the jet 
frame magnetic field probed by parsec scale VLBI observa- 
tions is dominated by either the azimuthal or axial field, 
then the transverse profiles described here would not ex- 
hibit the asymmetries that are the subject of this paper 
(except for the RM profiles). Ideal magnetohydrodynamic 
(MHD) models of jets launched from accretion disks predict 
that, in the conical region of the jet, the observer frame az- 
imuthal to axial field ratio at a distance z above the disk 
is B^/Bz « Rj{z)Q,/c, where is the angular frequency of 
the magnetic field line's footpoint anchored in the accretion 
disk at a radius vq from the central black hole, and Rj is 
the cylindrical radius of the jet at a distance z above the ac- 
cretion disk (e.g., Appl fc Camenzind||1993] ). Assuming the 
field line footpoint is located at ro = r}GM/(? and is in Ke- 
plerian orbit around a black hole of mass M, we estimate 
the following for the jet frame ratio for typical VLBI AGN 
jet parameters: 



B' 



(5) 



-3/2 



O.lpc 



f M 

V 109 Mq 



(1) 



This simple estimate for B'^/B'^ is sensitive to the jet pro- 
duction location, ro — parameterized by 77 in equation ([TJ — 
which has been estimated to be anywhere from r) ~ few to 
several tens depending on the spin of the black hole and 
other properties of the AGN system ( Meier et al.|2001 Vla- 
|hakis fc~K onigl 2004|. If this wide range of values for r] is 
realized in nature, then it is possible that a modest fraction 
of parsec and sub-parsec scale jets probed by VLBI exper- 
iments have magnetic fields where B'^/B'^ ~ 1, producing 
the asymmetric transverse profiles discussed in this paper. 

Alternatively, B'^/B'^ may not depend on conditions 
at the base of jets as it does in equation ([T]). If ultra- 
relativistic AGN jets are magnetically dominated on parsec 
scales ( |Blandford|2002[|Vlahakis fc Konigl|2004[ ) and condi- 
tions on the boundary of the jet do not change too rapidly, 
the jet may relax to a force-free structure, V x _B = kB 
( Choudhuri fc Konigl|1986 1. Force- free fields have the prop- 
erty that B^/B^ ~ 1, implying that significant transverse 
asymmetries are expected in VLBI profiles. Two simple 
force- free configurations to which a cylindrical jet could re- 
lax are (i) a diffuse pinch in which k depends on position 
(e.g., Lynden-Bell|1996" \ and (ii) the option used in this pa- 
per: a reverse field pinch where k = constant (Lundquist 



|T950l [Choudhuri fc Konigl|1986| ). The mmimum energy con- 
figuration of a helicity conserving magnetic field corresponds 
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to a force- free field in which k — constant (Woltier 19581. 
Thus, if the jet frame field relaxes via magnetic dissipation 
while approximately conserving helicity, the field may find 
the minimum energy state (Taylor 19741, which is the re- 



verse field pinch for a cylindrical geometry: 

B'{p,(t,,z) = Bo[0,Jiikp),Joikp)], 



(2) 



where Ji,2 are Bessel functions of the first kind, and the 
cylindrical coordinate, p, is normalized so p = 1 at the 
boundary of the jet. As can be seen in figure 1(b) we set 
k = 2.405, to ensure = at p = 1, such that there are 
no reversals of the axial field in the emission region, or jet 
spine. For the purposes of calculating the _RAf, we locate the 
Faraday rotating region outside the emission region between 
cylindrical radii p = 1 and p = 1.6. The decision to locate 
the Faraday rotating region outside the emission region is 
motivated by a variety of observations described in ^ 



2.2 Jet Geometry and Emission Properties 

We assume the jet is a steady cylindrical flow of radius Rj 
with a bulk Lorentz factor of F = (1 — /3^)~^''^ = 10, where 
P — v/c = /3(0, 0, 1) is the jet speed in units of the speed of 
light expressed in cylindrical coordinates (p, </!>, z) centered 
on the jet axis. The observer is located in the = (or 
x-z) plane such that observed photons move along the unit 
vector n = (sin^o;,, 0, cos^o^), where 9ob is the angle between 
the photon propagation vector, n, and the jet propagation 
direction, z. 

The cylindrical jet approximation used here breaks 
down when the jet is viewed with dob ^ dj, where 9j is the 
half opening angle of a more realistic conical jet. Thus, com- 
paring our model to blazars (AGN where Got l/F) might 
seem questionable except for the simple trend in VLBI ra- 



dio jet surveys such as the MOJAVE sample (Lister et al 
2009): a significant majority of MOJAVE sources have a 
core-jet structure suggesting that, typically, 9ob > 9j. A jet 
viewed with 6ob ^ 9j would not have a simple core-jet mor- 
phology. When viewed at such small angles, the jet would 
have features moving away from its core in all directions on 
the sky, leaving a core without a distinct jet morphology. As 
the VLBI images of most of the sources in MOJAVE reveal 
a clear core-jet morphology, we can conclude that they are 
viewed with Oob > dj- 

We assume the jet synchrotron emission is produced 
by a power-law distribution of relativistic electrons, dn' — 
K'i.E'~^dE' . The polarized emission coefficients in the 
plasma rest frame for a power-law distribution of electrons 
can be expressed as (leaving out the primes) 



j« =cWif.|Bsinx|"*'^-<''-^'/^ 
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Figure 1. (a) This schematic shows one magnetic field line in 
the jet comoving frame from the observer's perspective where 
e^jj < 7r/2 and B'^/B'^ ^ 1. On the right side of the jet sin x' ~ 1 . 
that is, the magnetic field is approximately perpendicular to the 
observer's line of sight. On the left side sinx' ~ Oi meaning that 
the magnetic field is pointing almost directly at the viewer. This 
asymmetry manifests itself in VLBI profiles of intensity, polar- 
ization, spectral index, and RM. Hereafter, the terms "left" and 
"right" will be referencing the sides of a jet viewed with 9'^^ < 7r/2 
that is filled with a right-handed helical magnetic field as shown 
here, (b) This is a plot of equation (j2]l as a function of cylindrical 
radius, p, which is the magnetic field configuration used in this 
work. Note that the pitch angle of the magnetic field configura- 
tion varies with cylindrical radius such that the field is purely 
axial on the jet axis and purely azimuthal on the jet boundary. 



and the absorption coefficients are 



(0 {») 



;^e|i3sinx|<''+'>/'j^-'^+''/' 
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(4) 



where x is the angle between the magnetic field and line of 
sight (i.e. Bcosx = B ■ n), is the electron mass, e is 
the elementary charge, v is the frequency of the electromag- 
netic wave, and T{x) is the gamma function of argument x 
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( Pacholczykl|1970[ ). The i = 1 (i = 2) state refers to the 
value of the coefHcients with the upper sign (lower sign) in 
equations (|3| and (j4|. These coefficients correspond to or- 
thogonal radiation polarization states for which the photon 
electric vector is perpendicular (parallel) to the component 
of the local magnetic field projected onto the sky (i.e. per- 
pendicular to the photon propagation direction). 

The asymmetry highlighted in this paper is that the left 
and right sides of the jet axis projected onto the sky have dif- 
ferent synchrotron emission properties (cf . |Aloy et al.|2000[ ) . 
The cause of this asymmetry is most clearly illustrated by 
the representation of a single magnetic field line viewed in 



the jet frame in figure 1(a) where the left side clearly ap- 



pears different from the right side of the helix. (Hereafter, 
when referencing the "left" or "right" side of the jet, we 
are referring to the emission region of a jet with a helical 
field that is a right-handed helix for which 9'^i, < n/2 as 
More precisely. 



on the 



shown in figure 1(a) ) More precisely, {\smx' l/i^ft 
left side of the jet is less than {\sinx'\) ^ight right 
side of the jet, where x is the angle between the jet frame 
line of sight, n' , and magnetic field, B' . The {) symbols re- 
fer to an emission weighted average along the line of sight 
intersecting the jet. As discussed in |2.1[ there is only a sig- 
nificant asymmetry in (jsinx'l) if B'^/B'^ ~ 1; otherwise, 
if the field is azimuthally dominated or axially dominated, 
then (|sinx'|)ie/t ~ il^^^x'l) right ^^"^ transverse asym- 
metry will be observed. This asymmetry in (|sinx'|) can 
clearly cause asymmetries in the jet's synchrotron properties 
as both the synchrotron emission and absorption coefficients 
depend on |sinx'| as seen in equations ([sjl and Q. 

Note that the asymmetry in {|sinx'|) is not intrin- 
sic, but instead depends on the jet frame viewing angle 
O'^i,, which is connected to F and 9ob via the relativistic 
aberration effect. For O'^i, — 7r/2 there is no asymmetry 
because {|sinx'|),^/t = (IsinxDrighf The calculated pro- 
files are asymmetric when the jet frame viewing angle is 
^o6 = ^/'^ i ^, where the constant ^ is between and 7r/2. 
Profiles originating from a jet with S^j, = 7r/2 -I- ^ or a jet 
with O'^i, = tt/2 — ^ are refiected versions of the other, where 
the reflective symmetry axis is the z-axis projected onto 
the sky. For example, if an AGN jet with ^^j, — 7r/2 + ^ 
and a positively skewed profile has its z-axis rotated so that 
^ob = ~ ^) then the rotated AGN jet's new profile will 
be negatively skewed (i.e. refiected about the z-axis pro- 
jected onto the sky). In a relativistic flow where F 3> 1, no 
asymmetry is expected for the observer frame viewing an- 
gle 6oh ~ 1/F, corresponding to O'^i, ~ n/2. Furthermore, 
the observer frame viewing angles Oot = N/T or l/{Nr), 
where N <ti T, correspond to O'^i, ~ 7r/2 ± ^. Consequently, 
AGN jets viewed with 9ob ~ N/T produce approximately 
the same profile as 6ob ~ 1/(A''F), except with the profile 
being refiected about the jet axis projected onto the sky. 



3 INTENSITY PROFILES 

In a large-scale magnetic field, optically thin synchrotron 
emission is not isotropic because it depends on sinx'. Since 
the synchrotron emissivity is an increasing function of sinx', 
more power is emitted toward observers whose line of sight 
is perpendicular to the magnetic field. Thus, for a helical 
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Figure 2. This plot depicts theoretical intensity curves calculated 
as described in section [s] for observer frame (jet frame) angles 
of 9,br = 1/1.2 (74°), 6»„6r = 1/2 (53°), e„i,r = 1/3 (37°), 
= 1 (170°), where T = 10. Except for the e^b = 1 curve, 
all of the curves are calculated using dob^ < 1, so that the jet 
frame viewing angle is 6'^^ < tt/2 and the jet is viewed as shown 
in figure |l(a)| Note that the intensity peaks occur on the right 
side of the jet (a > 0) where, as shown in figure |l(a)[ the angle 
between the line of sight and the B-field, x'l is closer to ~ it/2 
and therefore, according to equation the emission is greater. 
The exception to this trend is the O^b = 1 curve which represents 
jets that are close to being in the plane of the sky. Radio galaxies 
with relativistic flows behave differently because the rest frame 
viewing angle is almost 6'^^ ^ it. 



tensity will be observed on the right side of the jet (a > 0) 
where sinx' ~ 1 than on the left side where sinx' ~ 0. We 
demonstrate this behavior by integrating over the jet-frame 
emission coefficient defined in equation (|3| for a power law 
distribution of electrons: 



Ii, oc 



s' 



jl'ds', 



(5) 



magnetic field as seen in figure 1(a) more synchrotron in 



where S is the total path length through the jet for a given 
line of sight. 

The results of these calculations with arbitrary normal- 
ization can be seen in figure [2] Doppler beaming is not im- 
portant for the shape of these profiles as it only affects their 
normalization. As expected, the profiles go to zero at the 
edges of the jet and peak in the projected central region be- 
cause lines of sight are longer through the central regions 
of the jet than on the sides. The anisotropic synchrotron 
emission in the presence of the helical field is responsible for 
the off-center peaks and most importantly for the skewness 
of the profiles. Thus, the primary observable effect helical 
fields have on intensity profiles is to skew them such that 
the longer tail of the profile is on the side of the jet where 
the magnetic field is approximately parallel to n' . An im- 
portant transition occurs in the skewness of our profiles at 
dob ~ 1/F. Given the handedness of the helical field in fig- 
ure |l(a)[ a jet viewed with 9ob < 1/F has intensity profiles 
with a negative skew (i.e. its long tail is to the left), while 
the same jet viewed with 6ob > 1/F has a positive skew. 
The skewness is most clearly manifested for Bob < 1/2F and 
9ob > 2/F; otherwise the jet's intensity profile will be ap- 
proximately symmetric. 

As the 9ob ~ 1 case in figure [2] shows, the skewed profiles 
become double humped for either radio galaxies, 6ob ~ 1, 
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Figure 3. (a) This observed 5 GHz intensity profile of 3C 78, with 
a beam-size (FWHM) of 2.5 mas X 2.5 mas, displays none of the 
predicted asymmetry. This could be because predicted asymmet- 
rical features are below the sensitivity of VLBA. (b) This observed 
15 GHz intensity profile of the quasar NRAO 140 (with a beam- 
size of 0.9 mas X 0.6 mas) does exhibit a small asymmetry, a 
longer tail on the left, which is consistent with our theoretical 
profiles. 



or for blazars with 6ob ^ 1/r since in botli cases the rest 
frame viewing angles are 9ob ~ vr or respectively. This 
edge-briglitening occurs because for such jet viewing angles 
the jet frame line of sight is almost parallel to the jet axis. 
Therefore, lines of sight crossing through the center of the 
jet where the axial field dominates will be parallel to the 
magnetic field, and little synchrotron radiation will be emit- 
ted towards the observer. Closer to the edges of the jet the 
magnetic field is predominately toroidal and perpendicular 
to the line of sight, so more radiation is emitted toward 
the observer. Thus, the presence of edge-brightening in jets 



such as M87 ( Reid et al. 1989 1 and 3C 345 ( Unwin et al. 
1992[ ) may be explained by jet viewing angles with 6ob ^ or 

< 1/r. 

Observed intensity profiles of the radio galaxy 3C 78 



and blazar NRAO 14(Q are shown in figures |3(a)| and |3(b) 
respectively. 3C 78 is a radio galaxy with Got ~ 1 and a 
mildly relativistic parsec scale jet (/3 5C 0.52) as indicated 
by its parsec scale jet to counterjet surface brightness ratio 
( Kharb et a"L||2009| ). The symmetry of 30 78's profile may 
be explained by the fact that the jet frame line of sight is 
perpendicular to the jet axis, or it is possible that the finite 
beamwidth smears out any small skewness the true intensity 



profile contains. NRAO 140 (fig. 3(b) I is a relativistic jet 
exhibiting high pattern speeds of f3app ~ 13c and therefore 
9ob ^ 9° < |Lister et aL]|2009[ ). The intensity profile of NRAO 
140 displays an asymmetry in the form of a hump in one 
tail of the profile consistent with our theoretical intensity 
profiles. 

Further details regarding the observations of the inten- 
sity profiles and other observables covered in sections i|3j 
^ and ^ are as follows: The observations for 3C 78 were 
carried out with the VLBA and Effelsberg on September 10, 
2005 ( [Kharb et al.|2009[ ). NRAO 140 was observed as part of 
the MOJAVE program with the VLBA on August 9, 2007. 
The slices in total intensity (and fractional polarization, see 
figure[6| were obtained with the task SLICE in AIPS. These 
were obtained roughly perpendicular to the local jet direc- 
tion, and at a distance of approximately 5 mas for 3C 78, 
and 4.5 mas for NRAO 140. We note that the observed pro- 
files of total intensity (and fractional polarization) do change 
somewhat along the jet at different distances from the core. 
Therefore, to overcome this problem in future works, it will 
be necessary to study a statistically significant sample of jets 
and multiple slices within each jet. 



4 FARADAY ROTATION PROFILES 

A linearly polarized electromagnetic wave can be decom- 
posed into left-handed and right-handed circular polariza- 
tion states which propagate through a magnetized ionic 
plasma with different phase speeds, causing a rotation of 
the plane of polarization. This effect, known as Faraday ro- 
tation, can be calculated in the cold plasma approximation 
wherein the angle through which the EVPA is rotated is 



A- ^ \'2 



(6) 



where is the component of the jet frame magnetic field 
parallel to the jet frame photon propagation vector, n'j, is 
the jet frame number density of thermal electrons, and A' is 
the jet frame electromagnetic wavelength ( Burn|1966 1. 

The discovery of unambiguous RM gradients first in 3C 
273 ( |Asada et ar]|2002| [Zavala fc Taylor||2005|l and then in 



other AGN jets (e.g., [Kharb et al. 2009[ jCroke et aL][20T0 l 
is consistent with the prediction that a large-scale helical 
field will give rise to a smoothly changing RM across the jet 



(Laing 1981 Blandford 19931. Various parsec scale polari- 



metric observations suggest the Faraday screen responsible 
for RM gradients is not cospatial with the synchrotron emit- 
ting region (e.g., [Zavala fc Taylor[[2004[ [Kharb et al.|2009 |. 
The remaining probable locations of the Faraday screen are 
the sheath of the jet, the broad line region (BLR), and the 



1 NRAO 140 data obtained from the MOJAVE database at 



http: / /www.physics.purdue.edu/MO JAVE/allsources.shtml| 
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narrow line region (NLR). However, the BLR and NLR have 
been ruled out specifically in 3C 273 by the short time vari- 



ability of RM ( Asada et al.| 2008a) and in general by volume 
filling factor arguments and other polarimetric data ( Zavala 
&: Taylor|2004 \ . Thus, in this paper, we model the Faraday 
screen as a cylindrical shell of thermal particles (a sheath) 
that likewise carries the helical field defined in equation ([2|, 
but primarily resides outside of the synchrotron emitting 
spine of the jet. 

The radial thermal electron density profile function we 
have chosen is 

Wt oc p^exp (-p^), (7) 

within the emission region where p < 1 the thermal electron 
density goes as u't ~ and is therefore suppressed until 
p ~ 1 where it sharply rises. This ensures the Faraday rota- 
tion occurs primarily in a cylindrical shell surrounding the 
jet up to a cylindrical radius of p = 1.6 (see figure 1(b) I. 
Consequently, to match the observational evidence in our 
model there is little overlap between the emission region 
and the Faraday screen. We carry out the calculation for 
the Faraday rotation in the rest frame jet assuming that the 
Faraday rotating sheath is moving at the same speed as the 
jet (F = 10): 



Axf = a n'j- \B'\ cosx'ds' . 



(8) 



Figure 4(b) shows the calculated invariant quantity Axf ~ 
RMX^ from equation ^ for different viewing angles. The 
normalization, A, is not as important as the shape of the 
RM profile, but is set to give amplitudes of Axf more or 
less consistent with observations. As expected there is a gra- 
dient in RAI across the jet which is due to the line of sight 
component of toroidal magnetic field, B^, changing across 
the jet. The asymmetry in the amplitude of RM in the jet is 
also apparent: the \RM\ on the left side of the jet is less than 
the |_Ril/| on the right side. This is due to the axial field's 
contribution to the jet frame line of sight magnetic field. For 
viewing angles Oot < l/F the sheath magnetic field (see fig. 
1(b) I is aligned with the line of sight on the right side more 
than on the left side of the jet. 

The obser ved RM profile for 3C 78 ( [Kharb et al.|2009 1 



in figure 4(c) displays the asymmetry discussed above: the 
magnitude is greater on one side (the left) of the jet than it 
is on the other, suggesting the sheath magnetic field is more 
aligned with the line of sight on the left side of the jet. Other 
observed RAI profiles reveal behavior more similar to our 



calculated profiles. Asada et al. (2008b I, for example, detect 



an RM gradient using a different cut of NRAO 140 wherein 
the RM changes sign across the jet with the gradient of 
RM steepening toward the edges of the jet. Also of note, 



Asada et al. 



(2008b| finds that the magnitude of RM is 
roughly symmetric across the jet, suggesting that the axial 
field is not contributing to the RM. This symmetry in the 
magnitude of the RM suggests that either the jet frame line 
of sight is orthogonal to the jet axis (i.e. Oob = l/F) or the 
sheath field could be toroidally dominated. 



5 FRACTIONAL POLARIZATION PROFILES 

The linear fractional polarization, 11, primarily depends on 
the geometry of the magnetic field in the emitting region. 
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Figure 4. (a) The dotted line displays the density profile of non- 
thermal synchrotron emitting particles. The solid line shows equa- 
tion ^ which is the density profile of thermal particles in the 
Faraday screen. The small amount of overlap between the two 
profiles produces little internal Faraday rotation, (b) Theoretical 
profiles of Faraday rotation Axf ■ The primary effect viewing an- 
gle has on such profiles is to shift them up or down. The amplitude 
|i?Af I is higher on the right side of the jet because the magnetic 
field in the sheath (see fig. |l(b)| is more aligned with the line of 
sight on the right side of the jet than on the left, (c) This RM 
slice for 3C 78 was obtained where the RM gradient was clearly 
observed, at a distance of ~ 4 mas from the core. Note that the 
|/?A/| is higher on the left side. This implies the sheath magnetic 
field is more aligned with the line of sight on the left side than on 
the right side of the jet. The beam-size for this image is 2.5 mas 
X 2.5 mas; for further error analysis, see |Kharb et al.] ( |2009[ l. 



If the magnetic field is isotropic and disordered in regions 
smaller than the beam, then FT ~ 0, and if the magnetic 
field is uiuformly oriented, then the fractional polarization 
reaches a maximum of 11 = (p-|-l)/(p 4-7/3) ~ 0.7 for typical 
values of the electron index, p 



2-3 ( |Pacholczyk||l970 |. 
Thus, lines of sight that pass through the edges of a jet with 
a helical field will encounter a more uniform field than lines 
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of sight passing through the center of the jet, providing an 
explanation for the often observed rise in fractional polariza- 
tion towards the edge of jets (though it should be noted that 
observations of fractional polarization at the edges of jets of- 
ten have large errors). In observational transverse profiles, 
the linear polarization tends to rise towards the edges of the 
jet, and 11 is often significantly higher on one edge of the jet 



( Attridge et al.|1999||Pushkarev et aLpOOSlfZava la fc Taylor 
Gomez et al.||2008 |. This gradient in 11 across the jet 



2005 



(in which 11 is higher on one edge of the jet) can also be in- 
terpreted as a signature of a helical field (discussed below). 
Observed polarization profiles also tend to be more strongly 
asymmetric than intensity profiles as seen, for example, by 
comparing figure [3] to figure |6] below. 

To calculate the fractional linear polarization for a syn- 
chrotron emitting plasma with relativistic bulk motion we 



follow the procedure outlined in Lyutikov et al. ( 2003 2005 1 



and calculate the Stokes parameters I, Q, and U {V = 0): 

rS 

Qoc / (Bsinx')*''"^''^'cos2(x + AxF)ds 



)7(x f\Bamx'y^+''>^\m2{x + AxF)ds, (9) 
Jo 

where x is the angle the EVPA makes with the projection 
of the jet axis onto the sky (measured clockwise), and Axf 
is the angle through which the EVPA is rotated in the fore- 
ground Faraday screen (see [Q. If there were no bulk rel- 
ativistic motion involved, the observed x would be perpen- 
dicular to the projection of the jet magnetic field onto the 
sky. However, in our case where the emitting fluid element 
has a relativistic bulk velocity we incorporate the effects of 
relativistic aberration on x ( |Blandford fc Ko nigl|1979{ |Lyu-| 
tikov et al.][2003 , 2005]). Equations (|5|) and ([9| allow us to 
calculate the fractional polarization: 



n = 



(10) 



I 

We found the RM profile calculated in ^had no effect on 
the qualitative features of the profiles. For this reason we 
set Axf = in order to calculate the intrinsic EVPA angle 
which, given the cylindrical symmetry of the jet, can only be 



parallel or perpendicular to the jet axis on the sky ( Lyutikov 
et al.|2005| ). Thus, we plot 

Q 
I ' 



n 



(11) 



which displays both the linear fractional polarization {U is 
always zero due to cylindrical symmetry) and the EVPA 
direction. When 11 is negative, the EVPA is perpendicular 
to the jet because x = 7r/2, or Q oc cos2x < 0. When 11 is 
positive then the EVPA is parallel to the jet, so x = and 
Q oc cos2x > 0. 

The calculated fractional polarization profiles are shown 



in unconvolved form in figure 5(a) The magnitude of linear 
polarization reaches the maximum at the edges of the jet, 
decreases towards the center, then, after reaching zero, in- 
creases again in the center. (The increase in 11 in the center 
region would appear as a "bump" if the magnitude of 11 
were plotted.) This "bump" corresponds to an EVPA fiip 
from perpendicular to the jet a^cis on the edges to parallel 
to the jet in the middle region of the jet. These fiips have 
also been predicted by other helical models (e.g. Laing|1981 




-0.8 -0.6 -0.4 -0.2 0.2 0.4 0.6 0.8 
projected transverse distance, a 



(a) 




-1.5 -1 -0.5 0.5 1 1.5 
projected transverse distance, a 



(b) 



Figure 5. (a) Shown here are profiles of theoretical linear frac- 
tional polarization, 11 = Q/I, for difi'erent viewing angles. The 
sign of the polarization fraction refers to the EVPA direction: 
negative 11 corresponds to an EVPA that is perpendicular to the 
projected jet axis on the sky and positive 11 corresponds to an 
EVPA that is parallel to the projected jet axis on the sky. As ex- 
pected for large-scale magnetic fields, the polarization increases 
towards the edges of the jet. The most notable qualitative asym- 
metrical feature is that the left side of the polarization profile has 
a much steeper gradient than the right side, (b) This plot shows 
the same theoretical profiles convolved with a Gaussian beam. In 
these plots the primary effect is that the polarization of the left 
side of the jet is lower than the right side. 



Lyutikov et al.|2005|[Broderick fc McKinney|2010[ |. It is im- 
portant to note that this "bump" is not a robust feature 
of our calculations. That is, for other configurations of the 
magnetic field such as the diffuse pinch, such a "bump" does 
not appear for all viewing angles. The most robust feature of 
the calculated profiles is the difference in the gradient of lin- 
ear polarization across the jet. The polarization decreases to 
zero on the left side of the profiles much more quickly than 
on the right side. 



The profiles shown in figure 5(a) are unlikely to be di- 
rectly detected because the sharp features will be smoothed 
by a finite beamwidth. Thus we convolve each of the stokes 
parameters with a beam which has a standard deviation 
of a — 0.2 units (where the true projected jet diameter is 2 
units). We define the edges of the convolved polarization pro- 
files as the points at which the convolved intensity reaches 
0.5% of the intensity maximum. These points represent the 
lower limit of the detector. The results are shown in figure 
5(b) The most robust feature is a gradient in polarization: 
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Figure 6. The slices in fractional polarization (and total inten- 
sity, see figure |3| were obtained with the task SLICE in AIPS. 
These were obtained roughly perpendicular to the local jet direc- 
tion, and at a distance of approximately 5 mas for 3C 78, and 
4 mas for NRAO 140. The fractional polarization profiles also 
indicate the +/- sigma errors, (a) In this observed polarization 
profile of the parsec-scale jet of 3C 78, the large errors near the 
edge of the jet make it difficult to conclude much about the be- 
havior (such as the dip on the left), except that the polarization 
appears to higher on the right side of the jet (the beam-size is 2.5 
mas X 2.5 mas), (b) In this observed profile of NRAO 140 there 
is no prominent "bump" , but the fractional polarization is clearly 
higher on the right side. It also exhibits the expected increase in 
polarization towards the edges of the jet (the beam-size is 0.9 mas 
X 0.6 mas). 



the left side of tiie jet lias a lower polarization tlian tiie 
rigiit side. As in tlie unconvolved case, EVPA flips occur as 
tlie linear polarization passes tiirougli zero. Tiie number of 
flips in our calculated profile depend on the intensity level 
at which we cut off the polarization profile. (For example, if 
the n profile were cut off at 5% of the maximum intensity 
EVPA flips would be rarer.) Convolving the 11 proflles also 
moves the "bump" closer to the edge of the jet, and in the 
case of dob = l/l.Sr the "bump" is not apparent at all and 
the EVPA flips only once or not at all going across the jet. 



The most robust prediction from our convolved theoret- 
ical profiles, a gradient in 11 across the jet, is seen in 3C 78 
and NRAO 140 (fig. |6| and elsewhere in the literature ( At- 
tridge et alTjpgS), Pushkarev et al.|[2005l [Zavala fc Taylor 



2005 



Gomez et al.||2008[ ). We note that the observed pro- 



files of fractional polarization (and total intensity) do change 
slightly along the jet at different distances from the core. 
Also, just as in the observed profiles of intensity (fig.jsj and 
polarization (fig. [6|, the calculated unconvolved and con- 
volved polarization profiles display significantly more asym- 
metry than the (unconvolved) calculated intensity profiles. 
However, our model does not explain the significantly lower 
polarization level of 3C 78 (0% to 10%) compared to most 
of our theoretical curves and NRAO 140. This discrepancy 
can be due to a high degree of Faraday depolarization in the 
regions surrounding 3C 78's jet, or because of a significant 



disordered component to the jet magnetic fleld ( Cawthorne 
1993, Zavala fc Taylor 2004). 



et al.i 



EVPA flips are not seen in 3C 78 and NRAO 140. Two 
explanations for this absence are possible: (1) the magnetic 
field configuration in the emission region is different from 
that used in this work (e.g. a diffuse pinch field instead of 
a reverse pinch field), or (2) the EVPA flip occurs near the 
edge of the jet where the signal to noise ratio is too low for 
detection. VLBI EVPA flips as predicted in the unconvolved 
proflles are detected in some jets (Kharb et al., in prep). 



6 SPECTRAL INDEX PROFILES 

Typical AGN jet morphology consists of an unresolved 
bright core with a flat or inverted spectral index, a ~ — 1 to 
(Ii, oc t'""), and a steep spectrum jet with a spectral index 
of Q ~ 0.7 (e.g. |Zensus||1997| |. In some cases the transition 
from flat or inverted spectrum core to steep spectrum jet is 
thought to mark the transition between the compact opti- 
cally thick portion of the jet and the optically thin portion of 
the jet ( |Marsche"r||2009[ ) . Individual observed spectral index 
maps of AGN jets sometimes reveal more complicated be- 
havior such as asymmetries in the tranverse direction across 
the jet as seen in fig ure [t] (e.g., [Kharb et al.|2009|[O'Sullivan| 
[fc Gabuzda, 2009, Savolainen et al.|2008 1 . These gradients in 
the spectral index are usually thought to indicate that the 
jet is interacting with an inhomogeneous external medium 
(jet-cloud interactions). A large-scale helical magnetic field 
may also produce spectral index gradients. We explore two 
possible ways helical fields could affect the spectral index: 



anisotropic particle distribution functions (^6.1 1 and optical 
depth effects (jja2[). 



6.1 Anisotropic Distributions 

In many shock particle acceleration models the electron dis- 



tribution function is anisotropic (e.g., Lloyd & Petrosian 
[2000[ ). Particle in cell (PIC) simulations of shock acceler- 
ation occasionally find the electron index is dependent on 
the angle between the line of sight and the magnetic field, 
x' (Spitkovsky, private communication). We investigate the 
possibility that an anisotropy in the electron index p is pro- 
duced by the particle acceleration process or another un- 
known mechanism and is not isotropized too quickly by 
pitch-angle diffusion (ch. 12, Kulsrud 20051. If p changes 
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Figure 7. Observed spectral index map and profile for blazar 
0954+658 (figure from [O'SuUivan fc Gabuzda||2009] . Note the 
transition from a symmetric self-absorbed core to a region with a 
gradient as expected if the differential optical depth is inducing 
the spectral index gradient. The spectral index gradient is seen 
just upstream from where the jet bends, thus this gradient may 
also be caused by the jet interacting with its surrounding medium. 



with the pitch angle of the relativistic electrons (the pitch 
angle is the angle between the electron's velocity vector and 
the magnetic field), then the spectral index a will depend 
on x' ■ That is, for a given angle between the line of sight 
and the magnetic field, x', only the population of electrons 
whose pitch angle a — x' will be observed due to the rela- 
tivistic beaming of individual electrons. Therefore, the ob- 
served spectral index will be q = 1) /2 in an optically 
thin jet. 

The distribution we investigate changes linearly from 
one value (px) when x' = 7r/2 to another (py) when x' = 0: 



P(X') = ^+P||(1 

TT 

„ _ p(x') - 1 



TT 



(12) 



We obtain the spectral index profile in figure |8(c)] by assum- 
ing an optically thin jet and integrating equation ([5| with 
P ~^ p{x')- The spectral index is constructed by numerically 
evaluating equation ([5| for difi'erent ly, thereby constructing 
an optically-thin power-law spectrum. The range of values 
used for py and p± is the strong shock value of the electron 
index p — 2 ( Drury|Pl983 1 and a steeper spectrum of p = 4 
consonant with a radiatively cooled spectrum (p — >■ p + 1) 
with an injected electron index of Pinj ~ 3 seen in the cos- 
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Figure 8. (a) This figure uses a synchrotron spectrum from a 
homogeneous source with a power-law distribution to illustrate 
how optical depth effects can create a spectral index gradient. 
The dotted lines represent power fits to different regions of the 
spectrum. If the left side of the jet has a lower average optical 
depth than the right side of the jet then the left side of the jet will 
have o ~ 0.7 range while the right side of the jet will be in the a ^ 
—2.5 to range, (b) Spectral index profiles at different distances 
from the supermassive black hole. Note that such profiles have 
maximum changes of spectral index of Aomax 1. (c) Spectral 
index profiles from an anisotropic distribution. Note that such 
profiles have maximum changes of Aomax ~ 0.6 and are not very 
asymmetrical. 



mic rays. 



The results of these calculations are discussed in 



6.2 Optical Depth Effects 

We assume AGN jets consist of a compact optically thick 
core which gradually evolves into an optically thin jet. The 
existence of helical fields will modify this transition so that 



one side of the jet will become optically thin at a different 
core distance, z, than the other side of the jet. Therefore, 
there will be a region in the jet where one side has a different 
optical depth than the other side of the jet. The average syn- 



chrotron optical depth, r is defined as ( Rybicki & Lightman 



10 E. Clausen-Brown, M. Lyutikov, P. Kharb 



1979) 



[B smx ) ds. 



(13) 



0954+658 ( |Q'Sullivan fc Gabuzda(2009l ) which displays the 
expected behavior of a spectral index gradient caused by op- 
tical depth effects: the core is symmetrical and self-absorbed, 
and a gradient in spectral index is seen in the transition 
region between optically thick and optically thin. Alterna- 
tively, because the spectral index gradient occurs just up- 
stream from where the jet bends in 0954-1-658, the spectral 
index gradient may be the result of a jet-cloud interaction. A 
figuration as viewed'in figure 1 1(a) [the optical depth on the spectral ind ex gradient observed in 3C 273 (fig. l „Savolainen| 



Equation ( 13 1 reveals that the optical depth is an increas- 



ing function of sinx', implying that r is higher on the right 
side of the jet where sinx' ~ 1 than on the left side where 
sinx' ^ 1 (see fig. 1(b) I. Thus, for a magnetic field con- 



right side of the jet (a > 0) will always be greater than the 
optical depth on the left side of the jet (a < 0). 

To gain insight into how the synchrotron optical depth 
changes from optically thick to optically thin, we estimate 
how the following quantities scale with z: the electron den- 
sity, the magnetic field, and the length of the line sight pass- 
ing through the jet. While our jet is modeled as a cylinder, 
we derive the scalings as if the jet were conical. Conical 
expansion of the jet flow implies that the electron density 
scales as oc z~'^, and this expansion along with flux freezing 



implies the toroidal magnetic flelds scale as cx 



(If the 



magnetic field were primarily poloidal then B oc z~^; how- 
ever this difi^erence does not appreciably change our results.) 
The conical geometry implies that the line of sight distance 



scales as oc z. These scalings, along with equation 1 13 1 imply 
that 



J ccz 



(14) 



where k is the average absorption coefficient, j is the emis- 
sion function, and p {= 2.4) is the electron index. To calcu- 
late the spectrum we numerically solve the rest frame trans- 



fer equations (equation (3.66) of Pacholczyk 1970 1 for the 
Stokes parameters which include the polarized emission and 
absorption coefficients that depend on z as defined in equa- 



tion (14 1. After evaluating the radiative transfer equations 



and obtaining the total intensity numerically for different 
u, a spectrum is constructed over a one order of magnitude 
interval in frequency (e.g. 1 — 10 GHz). The spectral index is 
calculated by finding the best-fit linear slope of the log-log 
spectral interval using a reduced squares fit as illustrated in 
figure 8(a) 



6.3 Spectral Index Results 



The calculated profiles are shown in figures |8(b)| and |8(c)[ 
One notable difference between these figures is that the pro- 
file resulting from an anisotropic distribution function is not 
as asymmetric as that arising from optical depth effects, 
though this difference may disappear with a different form 



of equation ( 12 1. Two more general differences between op- 



tical depth effects and anisotropic distributions are: (1) op- 
tical depth effects produce a greater range in spectral index 
which include both inverted and steep spectral indices (the- 
oretically, a = —2.5 to ~ 0.7), as opposed to anisotropic 
distribution functions which only produce steep spectral in- 
dices (theoretically, a = (p — l)/2 = 0.5 to 1.5), and (2) 
anisotropic distribution functions can produce spectral in- 
dex gradients anywhere in the optically thin jet while optical 
depth effects give rise to such gradients only near the AGN 
core/optically thin jet boundary. 

We include an observed spectral index gradient in 



et al.|2008 1 is also consistent with a helical field inducing the 
gradient via optical depth effects. The spectra associated 
with various components of the parsec scale jet of 3C 273 
reveal a self-absorbed spectrum from the core to where the 
jet widens, at which point the south side of the jet (compo- 
nent B3) displays a self-absorbed spectrum while the north 
side (component B2) exhibits an optically thin spectrum. 

Unfortunately, observations of spectral index gradients 
alone are explained equally well by jet-cloud interactions. 
However, as will be discussed in ^ if such spectral in- 
dex gradients correlate with asymmetrical features of other 
VLBI observables, the case for helical field induced spectral 
index gradients is strengthened. 



7 DISCUSSION AND CONCLUSIONS 

We have shown that a cylindrical jet with an axially sym- 
metric large-scale helical magnetic field produces asymmet- 
rical transverse profiles for intensity, RM , fractional polar- 
ization, and spectral index. The profiles' asymmetries arise 
from changes in |sinx'| (or cosx' for RM) across the jet, 
where x' is the angle between the line of sight and magnetic 
field in the jet frame. Therefore, the degree of asymmetry in 
these profiles depends on the magnetic field structure (as- 
sumed to be force-free in this work), the jet viewing angle, 
and the bulk Lorentz factor. 

Unknown structural details of the jet-including bound- 
ary conditions, certain aspects of the functional form of the 
magnetic field, and the relativistic particle density profile- 
introduce significant uncertainty into all of our calculated 
transverse profiles. However, the qualitative features of the 
magnetic field used in this work (i.e. B^ = and Bz = max- 
imum on the jet axis and |B| decreasing with cylindrical ra- 



dius) are shared with a variety of analytic (e.g. Choudhuri 



& Konigl 1986 Lynden-Bell 1996 1 and numerical (e.g., |Li 
et al.||2006| |McKinney||2006[ ) models of AGN jets. To probe 
whether the asymmetrical features highlighted here are de- 
pendent on the specific structural form of the jet assumed in 
this paper, we have calculated transverse profiles (not shown 
here) using alternate boundary conditions, magnetic field 
structures, and density profiles. A force-free diffuse pinch 
field is tested, as well as several relativistic density profiles: 
n'^^i oc current''^ (ohmic dissipation), n^^j oc B'^ (equipar- 
tition), and a Gaussian density profile. These alternatives 
and their permutations affect specific qualities of the pro- 
files (discussed below), but they do not affect the robust 
asymmetrical features underlined in this work. 

Another significant unknown not treated in this work is 
whether a disordered field component would erase the pre- 
dicted asymmetries. Regarding polarization, for example, a 
disordered field component decreases the fractional polar- 
ization of a uniform field by the factor ~ Bq/{Bq + BJ), 
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where -Bp is the uniform field and Bd is the disordered field 
( Burn|l966 l. For intensity profiles, increasing B^/Bq gradu- 
ally erases the asymmetry because the disordered field com- 
ponent provides the same contribution to the synchrotron 
emission function on both sides of the jet. Consequently, as 
Bd/Bo is increased, the predicted asymmetries in this work 
would be diluted, though the exact ratio of B^/Bq at which 
the asymmetries become undetectable depends on how well 
the jet is resolved and the particular magnetic field structure 
of the jet. 

Although the aforementioned unknowns in AGN jet 
structure are significant, some conclusions can be made 
about the theoretical profiles calculated in this paper. For 
polarization profiles (Sj5|, the details of the jet structure un- 
fortunately dictate how many EVPA fiips occur across the 
jet and whether the polarization reaches zero. Ifowever, the 
significant gradient in polarization for all of the various jet 
structures discussed above (except for those with a signifi- 
cant disordered field component) is a robust feature of our 
theoretical polarization profiles. For intensity profiles (ij3| 
and spectral index profiles ([j6|, the qualitative features do 
not depend on the details of the jet structure. Intensity pro- 
files display a skewness that primarily depends on viewing 
angle and become bimodal for large or small viewing angles, 
0ob ^ or <C 1/r. These bimodal intensity profiles may ex- 
plain the observed edge-brightening in radio jets like M87 
(Reid et al.||1989). Regarding spectral index profiles, there 



are two mechanisms by which a helical field gives rise to 
spectral gradients: optical depth effects and anisotropic par- 
ticle distribution functions. Optical depth effects produce 
spectral index gradients that range from inverted to steep 
spectral indices only on the boundary between the AGN core 
and the optically-thin jet. Anisotropic particle distribution 
functions can induce gradients anywhere in the optically thin 
jet but have a more limited range of spectral indices, all 
of which are steep spectra. For Faraday RAI profiles ([|4|, 
many qualitative features do depend on specific, unknown 
features of the jet structure such as the presence of magnetic 
field reversals in the jet sheath and the functional form of 
the thermal particle density profile. The implications of this 
uncertainty in RM profiles is discussed below. 

As we have shown, a helical magnetic field can clearly 
produce emission patterns that are asymmetrical, however, 
other phenomena such as jet-cloud interactions or a curved 
jet can as well. Consequently, the ubiquity of asymmetrical 
profiles in parsec scale AGN jets alone lends little credence 
to the hypothesis that large-scale helical magnetic fields 
are present. An unambiguous signature of large-scale helical 
magnetic fields would be asymmetrical features in profiles of 
different VLBI observables that correlate with one another 
as specified in this work. The most straightforward corre- 
lation is between intensity profiles and linear polarization 
profiles: 

• Intensity profiles will be skewed such that the longer tail 
will be on the same side of the jet where the polarization is 
lower. 

Detecting skewed intensity profiles may be difficult since the 
predicted asymmetry will be significantly smoothed out by 
the beam. However, asymmetries in polarization profiles are 
easier to detect with a finite beam since the resulting profile 
will display a gradient in polarization. 



Comparing intensity and polarization profiles with spec- 
tral index profiles is not as straightforward. Ideally, all VLBI 
profiles should come from the same transverse cut on the jet, 
but care must be taken where spectral index gradients are 
present. The differing optical depths in a transverse cut with 
a spectral index gradient will induce asymmetrical features 
in polarization and intensity profiles that may be difficult to 
disentangle from the asymmetrical features induced by heli- 
cal fields. Thus, it may be necessary to examine polarization 
and intensity profiles in optically thin regions downstream 
of observed spectral index gradients. Assuming that (a) the 
same helical field is present at both the spectral index gra- 
dient cut as well as the intensity and polarization transverse 



cut, and (b) differing optical depth effects (S 6.2 1 are respon- 
sible for the spectral index gradient, then: 

• The side of the jet where the spectrum is more optically 
thin (or steeper) will be the side where the intensity profile 
has a long tail and where the polarization is lower. 

However, if the particle distribution function of the rela- 



tivistic electrons is anisotropic (^6.1 1, there will be a corre- 
lation between the direction of the spectral index gradient 
and other VLBI observables, but we cannot determine what 
kind of correlation should exist without a better understand- 
ing of particle acceleration mechanisms in parsec scale AGN 
jets. Correlations seen between spectral index gradients and 
other observables far away from the core in the optically 
thin region of the jet would support the anisotropic electron 
distribution function explanation. 

RAI profiles probe the surrounding jet sheath, not the 
jet spine, where synchrotron emission takes place. Unfortu- 
nately, the unknown nature of the magnetic field and ther- 
mal particles in the sheath make it difficult to identify robust 
features of RM profiles besides the existence of transverse 
gradients. Therefore, RM profiles may not correlate with 
profiles of intensity, polarization, or spectral index. How- 
ever, if the spine and sheath of the jet contain ascending 
and descending magnetic fiux as in magnetic tower models 
(e.g., Lynden-Bell| 199 6| there may be a correlation between 
observables that probe the spine and RM , which probes the 
sheath. Thus, comparing RM gradients with other VLBI 
observables may further illuminate a possible connection in 
the magnetic fields between the jet spine and sheath. 

In this work we have compared observed profiles for 
NRAO 140 and 3C 78 with our theoretical predictions. For 
NRAO 140 the observe d inte nsity profile (fig. |3(b)[ ) and the 
polarization profile (fig. 6(b) I show the correlation expected 
from helical fields: the long tail of the intensity profile is on 
the side of the jet (the left) where the polarization fraction is 



low. 3C 78 is more ambiguous. Its intensity profile (fig. 3(a) I 



is not significantly skewed, but the polarization profile (fig. 
|6(a)[ ) does exhibit a gradient with higher polarization on the 
right side of the jet, though even this is unclear due to the 
large errors on the jet's edges. The spectral index map of 3C 
78 (not reproduced here, see fig. 4 of [Kharb et al.|[2009[ ) is 
difficult to interpret as it reveals gradients in both directions 
in different parts of the jet well beyond the optically thick 
core. If the polarization profile of 3C 78 is due to helical 
fields, a comparison can be made with its RM profile shown 
in figure 4(c) The left side of the 3C 78's jet exhibits a 



higher |-RM|, suggesting that 3C 78's magnetic field in the 
spine (i.e. emission region) and sheath is different from the 
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assumed magnetic field structure assumed in this work (fig. 
1(b) I. That is, the lines of sight through the left side of 3C 



78's jet spine and sheath are more closely aligned with the 
magnetic field than those on the right side. 

To test whether the signature of large-scale helical fields 
exist, a systematic search for the correlations between differ- 
ent VLBI profiles described here is required. This not only 
involves investigating as many resolved AGN jets as possi- 
ble, but also involves examining multiple profiles at different 
distances from the core for a single jet. In the future we in- 
tend to carry out such a search in the MOJAVE sample of 
radio loud jets with a focus on well-resolved jets. 
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